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Abstract 

The fidelity of laser-driven quantum logic operations on trapped ion qubits tend to be lower than 
microwave-driven logic operations due to the difficulty of stabilizing the driving fields at the ion 
location. Through stabilization of the driving optical fields and use of composite pulse sequences, 
we demonstrate high fidelity single-qubit gates for the hyperfine qubit of a m Yb + ion trapped 
in a microfabricated surface electrode ion trap. Gate error is characterized using a randomized 
benchmarking protocol, and an average error per randomized Clifford group gate of 3.6(3) x 10 -4 
is measured. We also report experimental realization of palindromic pulse sequences that scale 
efficiently in sequence length. 

PACS numbers: 03.67.-a, 03.67.Ac 


1 



The trapped atomic ion qubits feature desirable properties for use in a quantum computer 
such as long coherence times |TJ , high fidelity qubit measurement [2], and universal logic 
gates [3]. The quality of quantum logic gate operations on trapped ion qubits has been 
limited by the stability of the control fields at the ion location used to implement the gate 
operations. For this reason, the logic gates utilizing microwave fields HH6] have shown gate 
fidelities several orders of magnitude better than those using laser fields [7H9] . The UV laser 
beams used to drive Raman gates for a hyperfine ion qubit present a major challenge as they 
are subject to severe wavefront distortion in air due to turbulence, leading to amplitude and 
phase fluctuations at the ion location. 

Microfabricated surface electrode ion traps, where atomic ions are trapped above a two 
dimensional surface of electrodes, can provide a scalable platform on which to build an 
ion-based quantum computer [10, [IT]. Experiments using surface traps have demonstrated 
coherence times of more than 1 second [12], state detection with fidelities greater than 
99.9% [2], and low error single-qubit gates (e < 2.0(2) x 10 -5 ) using integrated microwave 
waveguides HE]. Use of high power UV lasers close to the trap surface can lead to substantial 
charging due to unwanted exposure pj. The recent development of single-mode fibers 
capable of delivering high power UV laser beams ra opens the possibility of significantly 
reducing the free-space UV beam path length and delivering a clean spatial mode to the 
ions, eliminating unwanted scattering off nearby trap structures. 

Here, we demonstrate low-error single-qubit gates performed using stimulated Raman 
transitions on an ion qubit trapped in a microfabricated chip trap. Gate errors are mea¬ 
sured using a randomized benchmarking protocol pa usi ns], where amplitude error in the 
control beam is compensated using various pulse sequence techniques [17, |18]. Using B2 
compensation HZ], we demonstrate single-qubit gates with an average error per random¬ 
ized Clifford group gate of 3.6(3) x 1CT 4 . We also show that compact palindromic pulse 
compensation sequences (PDn) [T8] compensate for amplitude errors as designed. 

Two hyperfine ground states of the lfl Yb + ion (|0) = 2 Si/ 2 |E = 0, m/ = 0) and |1) = 
2 S 1 / 2 1-f -1 = l,m/ = 0), shown in Fig. [IJi, separated by f qu ut ~ 12.6 GHz, serve as our qubit 
states. The energy separation between these states is relatively insensitive to the magnetic 
field fluctuations at the relevant magnetic field (~ 3 kHz/Gauss at 5 Gauss). Continuous 
wave (CW) lasers are used to perform Doppler cooling, resonant scattering for qubit state 
detection, and optical pumping out of the 2 D 3 / 2 state (not pictured in Fig. [1^). For Ra- 
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man transitions, we use picosecond pulses from a mode-locked titanium-sapphire (Ti:Sapph) 
laser frequency-doubled to a center wavelength of 376 nm, which creates combs in the fre¬ 
quency domain with comb teeth spacing equal to the laser repetition rate (/ rep ~ 76 MHz). 
The frequency doubler output is split into two nearly co-propagating frequency combs us¬ 
ing a single acousto-optic modulator (AOM2) driven with modulation frequencies f\ and 
/ 2 , as shown in the inset Fig. [D>- Resonant transitions are driven by pairs of optical fre¬ 
quency comb teeth ((2) and (3) in Fig. [I]d), one from each comb, with a frequency difference 
equal to f qu ut [121 El- As the repetition rate / rep of the laser drifts, the frequency differ¬ 
ence between these two comb teeth is actively stabilized to match / qubit by adjusting the 
modulation frequency /2 to maintain the quantity 166 x / rep + / 2 constant [20112T] . After 
the AOMs, the Raman beams are delivered through a 3-m single-mode fiber to a micro¬ 
electromechanical systems (MEMS)-based laser beam steering system that is capable of fast 
(< 2/xs), re-configurable addressing of individual ions in a linear chain with low crosstalk 
(< 3 x 10~ 4 ) [22J. The Raman laser power delivered to the MEMS system is actively 
stabilized using a gated digital proportional-integral (PI) loop (see below) [20] , After the 
addressing system, the co-propagating linear polarization of the Raman beams is converted 
to circular polarization with a quarter-wave plate to drive a cr + -transition. A dichroic filter, 
which reflects 376 nm light and transmits 370 nm light, is used to fold the Raman beams 
into a high numerical aperture (NA = 0.6, PhotonGear) lens that serves the dual purpose of 
imaging the ion fluorescence and focusing the frequency combs onto the ion. For state de¬ 
tection, the ion is imaged at 200 x magnification onto a photomultiplier tube (PMT) array. 
All experimental timing and measurement recording is carried out by a field-programmable 
gate array (FPGA) located in our main controller (shown in Fig. [Tj). 

For each experiment, the ion is first Doppler cooled for 1 ms using light that is red- 
detuned from 2 Si/ 2 | F — 1)—)- 2 Pi/2 |F = 0) resonance. The qubit is then initialized to the |0) 
state by applying CW light resonant with 2 Si/ 2 \F — 1) —)- 2 Pi /2 |E = 1) transition for 20 /xs. 
Following initialization, the intensity-stabilized Raman beams are pulsed on and off using 
the AOMs for a duration corresponding to the gate implementation. To measure the qubit 
state, light resonant with 2 Si/ 2 |E = 1) —> 2 Pi/ 2 |A = 0) transition is turned on for 400 /xs 
while ion fluorescence is measured using the relevant pixel of the PMT array. The ion will 
fluoresce if the qubit is in the |1) state, while the |0) state remains dark. 

Raman beam intensity fluctuations at the ion location result in systematic amplitude 
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errors which decrease gate fidelity. To minimize intensity fluctuations at the ion, a power 
stabilization system is implemented [2D], where a small amount of collimated light is picked- 
off after the beam exits the fiber and is incident on a photodiode (PD). The PD provides 
the error signal for a proportional-integral (PI) feedback loop that controls the amplitude 
of the direct digital synthesizer (DDS) signal driving AOM 1, as shown in Fig. [I] The PI 
loop is gated by a digital trigger pulse sent from the FPGA in our main controller. The 
Raman beams and PI loop are turned on during the Doppler cooling cycle, allowing the 
power to be stabilized. The loop is then turned off for the remainder of the experiment. 
Using this approach, the intensity fluctuations due to beam pointing instability before the 
fiber are corrected. Beam pointing instability after the fiber, due to air currents, is passively 
suppressed by enclosing the entire experiment in a box. 

A composite pulse sequence can be used in place of a single Raman pulse to make it robust 
against systematic errors such as amplitude, timing, crosstalk, or detuning errors nn HU [231- 
El and time-dependent control errors [2DJ. In our experiment, the impact of residual system¬ 
atic amplitude errors in the Raman beams is suppressed through the use of compensating 
pulse sequences. Since these techniques are usually designed to work on systematic errors 
that are constant over the duration of the sequence, the sequence length determines the 
bandwidth below which the effect of fluctuating error is suppressed citeKabytayev2014. The 
length of most compensating pulse sequences increases rapidly at higher error correction or¬ 
der. The palindromic pulse sequences (PDn) are unique in that they scale linearly with the 
corrected error order (to n = 12) [18]. Here we analyze the use of B2 [T7] and PD6 [IS] com¬ 
posite pulses and their ability to correct static amplitude errors in the presence of additional 
phase and timing errors for Clifford group gates. 

In the absence of noise, the target rotation R(9, 0) rotates the Bloch vector by an¬ 
gle 9 around the axis = X cos 0 + Y sin 0, represented by a unitary propagator 
R(9, 0) = exp (— \0a^). The B2 compensation sequence (also known as BB1), introduced by 
Wimperis na, is designed to correct the errors in the pulse area (due to amplitude or timing 
errors) to 0(e 2 ), where e is the fractional error in the control signal [23]. B2 compensation 
translates each single rotation into a sequence of 4 rotations. A target rotation R(9 t , 00 
becomes 

R(6 t , 00 R(9 t , 00-R(7t, 0t + 0 _b2)R(2vt, 0 t + 30 B2 )R(vr, 0 t + 0 B2 ), (1) 
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where 

^2 = cos- 1 . (2) 

A single B2 compensated pulse requires a total rotation angle of 9 tota i = 47r + 6 t , which 
requires time t to t a i = where Q is the Rabi frequency. While the B2 sequence is fairly 
short at n — 2, the Bn sequences increase in length exponentially with increasing n, requiring 
0(exp(n 2 )) pulses to reduce amplitude errors to 0(e n ) [2Tj . 

In comparison to Bn sequences, palindromic compensation sequences (PDn) scale effi¬ 
ciently in length, requiring only 2 n pulses ( 9 tota i = 2nn + 9 t ) to cancel errors to n th order, 
up to n = 12 [T8j. Here we use PD6 (n — 6), where a target rotation R(9 t , (fi t ) is replaced by 

R(9 t , (fit) R(9 t , <fit)R{ 7T, (fit + (fipDG-.l) R{n, (fit + 0PD6:2)R(tT, (fit + (fipDfr.3) R{n, (fit + (fipD6:4) 
R(n, (fit + (fiPDQ-.b) R{j( •) (fit + (fipD6:6)R(n : 4>t + (fipD&ti)R{n , (fit + <fi P DG:b) R{^ 1 + 4 >PDQ:a) 

-R( 7r ) & + (fipDG:3)R{T (j (fit. + (fiPD 6 . 2 )R(,7( 1 (fit + (fipD6:l) (3) 

with (fippe-.k given for all fc’s for <fi t = tt and 7 t/ 2 in Table [TJ 

The direct impact of Raman laser intensity fluctuation is to modify the amplitude of 
the qubit rotation, which is calibrated prior to each data set by fitting 3-5 periods of Rabi 
flopping. Raman laser intensity fluctuations also produce detuning errors due to a differential 
Stark shift between the two qubit levels. To limit the differential AC Stark shift, the power 
in each Raman beam is made roughly equivalent on the photodiode (PD in Fig. [Tj) by 
changing the amplitude of f\ and / 2 to ensure equal fiber coupling of both Raman beams. 
The hyperfine frequency, modified by the differential Stark shift, is found using a Ramsey 
interference experiment. We begin by initializing the qubit to the |0) state. The qubit is 
then placed into a superposition state using a 9 — 7t/ 2 microwave Ramsey pulse. This is 
followed by a wait time of up to 20 ms, during which Raman beams, 1.5 MHz off-resonant 
from the qubit frequency, are turned on allowing the differential Stark shift to modify the 
qubit frequency. Then the qubit is rotated again with a second 9 — tt/2 microwave Ramsey 
pulse. The final state of the ion will oscillate with the wait time when the microwave is 
off-resonant from the qubit frequency. The exact hyperfine qubit frequency in the presence 
of the Raman beams is found by scanning the microwave frequency to where the ion state 
is rotated completely to |1) at all wait times. This microwave frequency is used to compute 
the driving frequency f\ for AOM 2. Although we carefully calibrate the differential Stark 
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shift, drifts in the individual Raman beam amplitudes cause small changes in the qubit 
frequency (< 100 Hz), which is not effectively compensated by the pulse sequences used in 
our experiments. 

In a randomized benchmarking experiment, the qubit is first initialized to the |0) state, 
followed by a sequence of L gate operations, chosen uniformly and randomly from the 24 
Clifford gates (shown in Table [T|. A final Clifford gate is then chosen to bring the resulting 
qubit state to either the |0) or |1) state, at random. Then, the qubit state is measured and 
compared to the expected state. For each sequence length L , 20 random sequences were 
created and each sequence was measured 800 times, and the fraction of events where the 
measured result matched the expected result was recorded as the survival probability. 

The averaged survival probability per gate length is fit to the zero-order decay model [16] 

Fseq(L) = Aop L + Bq , (4) 

where F seq (L) is the survival probability at length L, p is related to the average error 
per Clifford gate (average error = ^), and A 0 and B 0 contain the state preparation and 
measurement (SPAM) errors and the error on the final rotation. 

The error on our measurement was calculated to account for both the variance due to 
the projective measurement statistics and the variance arising from the spread in fidelities 
of the underlying distribution of gate sequences, as outlined in [15] . An initial, unweighted, 
non-linear least squares fit with estimated SPAM errors (Ao = 0.47 and Bq = 0.517) was 
used to gain a first estimate of the underlying error, and from this an upper bound on 
the variance is calculated. This estimated variance is used to weight a second non-linear 
least-squares fit with floating SPAM parameters, and the resulting co-variance matrix was 
used to calculate the uncertainty. Figure [2] shows the result of a randomized benchmarking 
experiment for a single-qubit gate using the B2 compensating pulse sequence. Black dots 
represent the average survival probability of each random sequence measured. The average 
survival probability of all randomized sequences for each sequence length L (blue squares), 
with an upper-bound variance calculated as in Ref. [[15], are fit to equation (|4]) (red line) to 
extract an average error per Clifford group gate of 3.6(3) x 10 -4 . The error on the fit (light 
red) accounts for the distribution of measurements and the number of sequences measured 
for each L [15] . 

While we expect the actual detuning of the resonant combs arising from the change in 
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the differential AC Stark shift due to amplitude drifts of the Raman beams to be small 
(< 100 Hz), a much larger effective detuning error exists. The co-propagating Raman 
beams are circularly polarized, thus, other pairs of comb teeth (®-(2) pair, (3)-@ pair, 
and (T)-(4) pair in Fig. [1J>) can also induce Raman transitions, detuned by 5' ~ 4.5 MHz 
(for the case of ®-(2) pair, and (3)-@ pair) and 9 MHz (for the case of pair) in our 

setup. Further off-resonant beat-notes are also present, but have small contribution to gate 
errors. The contribution of these off-resonant Raman transitions add additional rotation 
to the desired state evolution, which is calculated by considering additional terms in the 
interaction Hamiltonian that describes the qubit subject to the driving held 

Hj = —{(Af cos 0 + Y sin 0) + 2[X cos (5't + 0i) + Y sin (S't + 00] 

+ [Afcos(2<W + 02) + Y sm(28't + 0 2 )]}, (5) 

where the rotations caused by the detuned Raman transitions due to the ®-(2) and (3)-@ 
pair act coherently on the qubit, and 0i and 0 2 describe the relative phases between the 
detuned Raman transitions and the resonant transition, which fluctuate with drifts in the 
repetition rate of the laser. The unitary propagator that describes time evolution of the 
qubit subject to this time-dependent interaction Hamiltonian can be computed using the 
standard Magnus expansion technique [28]. We use the second-order Magnus expansion to 
compute the time evolution operator for the qubit subject to each optical pulse. Once the 
fast oscillating terms (at 4.5MHz and 9MHz) are averaged out, the net impact of the off- 
resonant terms result in an effective time evolution operator R(9, 0) = exp (— \0\(T^ + 5Z]), 
where 5Z represents a deviation of the qubit rotation axis from the X — Y plane. 

For comparison with the experimental data, we simulated the exact benchmarking se¬ 
quences used in the experiments in the presence of these imperfections. To simulate each 
sequence of gates, we calculate the imperfect propagator R s (9,<f>) that describes the state 
evolution of each pulse, where 9 represents the rotation angle and 0 represents the phase 
of the driving pulse. Increase in pulse area due to static amplitude noise or timing error 
in the Raman beams leads to an error in the rotation angle of 6 — y 9(1 + e), where the 
multiplicative parameter e is the same for each pulse in a composite pulse sequence. The 
presence of off-resonant Raman beams and the (smaller) effect of the differential AC Stark 
shift are captured in the qubit detuning error S. To simulate a wider range of errors in the 
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pulse area, we purposefully change the timing of the pulses yielding the propagator 

R{0,4>) =exp[-^0(l + e)(<7* + <$Z)]. (6) 

Our calculation shows that 6 is effectively a random variable in the range of 0-3 kHz due 
to the drift in the repetition rate of the laser. We calculate a single propagator for each 
randomized benchmarking sequence by matrix multiplication of the imperfect propagators 
representing each individual pulse. The resulting final propagators are used to compute the 
final Bloch state, which is then compared to the expected state yielding a simulated survival 
probability. The simulated survival probability is fit to equation Q to produce an average 
error per gate for each series of sequences given specific noise parameters, which is then 
compared with the experimental results. 

Figure [3] shows a comparison of the error from uncompensated (primitive) gates and 
gates compensated using B2 and PD6 pulse sequences as a function of timing error. The 
B2 sequences can keep the gate errors to below 1% for timing error values in the range 
of | e |< 0.4, while the PD6 sequences can maintain similar gate error levels for timing 
error values of up to | e |< 0.6. Our simulations quantitatively reproduce the range over 
which the timing errors are compensated for both B2 and PD6 sequences using the imperfect 
propagator computed from the interaction Hamiltonian given in equation (J5]) . The minimum 
error value per gate is achieved by B2 sequence in our setup, since the longer sequence length 
of PD6 makes it more susceptible to the additional errors contributed by off-resonant Raman 
transitions. 

By driving the gates using two beams with polarizations that are orthogonal to each 
other and the quantization axis, the intra-comb contributions for the off-resonant Raman 
transitions (pairs ®-(2) and (3)-@) can be eliminated. The remaining leading-order error 
will come from inter-comb beat notes (pair ®-@) that can be further detuned by adequate 
choice of the repetition rate of the mode-locked laser. Our simulations suggest that these 
modifications can reduce the average error in the single-qubit gates by over one order of 
magnitude, at which point the gate performance will be limited by other detuning errors 
(such as the differential AC Stark shift due to amplitude drifts). 

In this work, we report high fidelity single-qubit gates driven with tightly focused laser 
beams on trapped ion qubits by laser intensity stabilization and use of compensating pulse 
sequences. An error probability as low as 3.6(3) x 10 -4 is demonstrated [29], consistent 


with error levels required for realizing a range of quantum error-correction schemes [301 - 132] . 
We experimentally verified the value of novel, length-efficient pulse sequences (PDn) that 
suppress errors to higher orders with modest sequence lengths. 
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4>t 


4>PD6:2 

</lPD6:3 

<t>PD6A 

<APD6:5 

4>PD6:6 

7 r 

0.38266 

-2.51430 

-1.75192 

0.05941 

2.67572 

0.39344 

7T / 2 

0.34769 

-3.06979 

1.55852 

-0.70890 

3.09692 

-0.62174 


TABLE I: PD6 Phase Angles 
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Clifford gate 

Physical gates 

Clifford gate 

Physical gates 

1 

I 

13 

Z/2 & X 

2 

X 

14 

X & Z/2 

3 

Y 

15 

Z/2 & X/2 

4 

Z 

16 

Y/2 & Z/2 

5 

X/2 

17 

X/2 & -Z/2 

6 

Y/2 

18 

Y/2 & Z 

7 

Z/2 

19 

-X/2 & Z/2 

8 

-X/2 

20 

-Z/2 & Y/2 

9 

-Y/2 

21 

Z & Y/2 

10 

-Z/2 

22 

-Z/2 & X/2 

11 

Z & X/2 

23 

X/2 & Z/2 

12 

X/2 & Z 

24 

-Y/2 & -Z/2 


TABLE I: Clifford group gates written as the physical gates (pulses) applied. 
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FIG. 1: Experimental setup. Inset (a) shows the relevant energy levels used in the 171 Yb + 
ion. Raman beams from a frequency-doubled picosecond Ti-Sapph laser passes through AOM 2 
driven by two frequencies (/i and / 2 ), and creates the two optical frequency combs as shown in 
inset (b). See text for a detailed description of the setup. DDS: Direct digital synthesizer, ADC: 
Analog-to-digital converter, PD: Photodiode. 
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FIG. 2: Randomized benchmarking of a single-qubit gate using B2 compensating pulse sequence, 
with average error per Clifford group gate of 3.6(3) x 10 -4 . 
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FIG. 3: Single-qubit gate error degradation with systematic amplitude error for primitive (black 
diamonds), B2 compensated (red squares), and PD6 compensated gates (blue circles). For some 
amplitude error values, multiple series of randomized sequences were measured resulting in multiple 
simulation lines over some amplitude error values. 
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